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ABSTRACT 

The hydro dynamic formation of massive black holes (BHs) in globular clusters is considered. In par- 
ticular, we examine the possibility of the BH formation induced by the radiation drag that is exerted 
on the interstellar matter by stellar radiation in globular clusters. The radiation drag extracts angular 
momentum from interstellar gas and thus allows the gas to accrete on to the cluster center. By incorpo- 
rating the realistic chemical evolution of globular clusters, we scrutinize the efficiency of the radiation 
drag to assess the total accreted mass. As a result, we find that if a globular cluster is more massive 
than 6 X lO^M© in the present-day stellar component, a massive BH with > 26OM0 can form within 
it. But, the BH-to-cluster mass ratio is considerably smaller than the BH-to-bulge mass ratio (~ 10~^) 
found in galactic bulges. The results are not sensitive to the assumed stellar initial mass function and 
star formation rate in the cluster, as long as the resultant color-magnitude relation, metallicity, and 
mass-to-luminosity ratio satisfy those observed in globular clusters. Hence, the putative linear relation 
between BH mass and bulge mass (Mbh — ^buige relation) cannot be extrapolated to globular cluster 
systems. In the present regime, we discuss the BH formation in M15, Gl, co Cen, the M33 nucleus, and 
the compact X-ray sources in M82. Finally, we argue observational indications of the formation process 
of massive BHs in globular clusters. We find that the final phase of BH growth due to the radiation drag 
can be observed as ultraluminous X-ray sources (ULXs) with ~ lO^^ergs"^. 

Subject headings: black hole physics - globular clusters: general - hydrodynamics - radiation 
mechanisms: general 



1. INTRODUCTION 

The recent compilation of the kinematical data on galac- 
tic centers has revealed that a central "massive dark ob- 
ject " (MDO), which is a candidate for suppermassive black 
hole (BH), does correlate with the mass of a galactic bulge; 
the BH-to-bulge mass ratio is ~ 0.002 as a median value 
(e.g., Kormendy & Richstone 1995). This correlation sug- 
gests that the formation of a supermassive BH is physically 
connected with the formation of a galactic bulge. 

If the BH-to-bulge relation can be extrapolated to small 
spheroidal systems like globular clusters, massive BHs with 
10^~^Mq may inhabit globular clusters. To date, some 
candidates of massive BHs in globular clusters have been 
reported; the estimated mass of a possible BH in M15 
is Mbh = (l-^t^.r) x lO^M© (Gerssen et al. 2003), or 
Mbh ~ 2 X lO'^M© in Gl (Gebhardt, Rich & Ho 2002). 
On the contrary, recent sophisticated numerical simula- 
tions (Baumgardt et al. 2003a, b) show that a massive 
BH with > 500 — lOOOM© is not mandatory to account 
for the observational data in Ml 5, and the luminosity pro- 
file of Gl can be well fitted without a massive BH. On the 
other hand, M33 nucleus does not appear to possess a mas- 
sive BH with an upper limit of ~ lO^M© (Gebhardt et al. 
2001; Merritt, Ferrarese & Joseph 2001). In addition, lat- 
est radio observations have suggested that the BH mass in 
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CO Cen should be less than about lOOM© for the spherical 
Bondi-Hoyle accretion rate, but the data are marginally 
consistent with a BH of about lOOOM© for more plausible 
accretion rate (Maccarone et al. 2005). Thus, the exis- 
tence of massive BHs in globular clusters is under debate. 

Some theoretical models have been proposed for the for- 
mation of massive BHs. One of them is the merger of 
multiple small BHs or stars (e.g., Lee 1987; Quinlan & 
Shapiro 1990). However, this scenario predicts that, even 
in a very compact star cluster like a globular cluster, the 
core-collapse is halted by the binary heating before the 
stellar density becomes high enough for stars to merge in 
a runaway fashion (Hut et al. 1992). However, it is pointed 
out recently that the mass difference of constituents plays 
an important role. Stars of different masses are not always 
able to reach the energy equipartition. This causes the 
heaviest stars to undergo core collapse on a time scale that 
is much shorter than the core collapse time for the cluster 
as a whole. Portegies Zwart & McMillan (2002) show by 
A^-body simulations that a runaway merger among these 
massive stars leads to the formation of a massive BH. On 
the other hand. Miller & Hamilton (2002) examine a model 
where a single BH (> 5OM0) located at the center of the 
cluster grows in mass through merging with stellar mass 
BHs. 

Another possibility is the hydrodynamic formation of 
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massive BHs through super massive stars. Recently, as a 
potential hydro dynamic mechanism in a bulge, Umemura 
(2001) has considered the effects of radiation drag, which 
is equivalent to a well-known Poynting- Robertson effect in 
the solar system. The radiation drag extracts angular mo- 
mentum from interstellar gas and thus allows the gas to 
accrete on to the center. This is a parallel model of the for- 
mation of massive BHs by the Compton drag in the early 
universe (Umemura, Loeb, & Turner 1997). The mass ac- 
cretion by the radiation drag has been explored in detail 
(Umemura, Fukue, & Mineshige 1997; Fukue, Umemura, 
& Mineshige 1997). The angular momentum loss rate 
by the radiation drag is given by dlnJ/dt —xE/c^ 
where J is the total angular momentum of gaseous com- 
ponent, E is the radiation energy density, and x is the 
mass extinction coefficient. Therefore, in an optically-thin 
regime, dhiJ/dt :^ —{tL^/c^M^q^^)^ where r is the to- 
tal optical depth of the system, is the total luminos- 
ity of the spheroidal system and Mgas is the total mass 
of gas. In an optically thick regime, the radiation drag 
efficiency is saturated due to the conservation of photon 
number (Tsuribe & Umemura 1997). Thus, an expression 
of the angular momentum loss rate which includes both 
regimes is given by dinj/dt :^ —{L:^/c^Mg){l — e~^). 
Then, the mass accretion rate is estimated to be M = 
—MgdlnJ/dt = L*/c^(l — e~^). In an optically-thick 
regime, this gives simply M = L^jc? (Umemura 2001). 
Then, the total accreted mass on to the MDO, Mmdo, is 
maximally Mmdo — j L^/c^dt. Although this is a simple 
estimation in one-zone interstellar medium (ISM), the ISM 
is observed to be highly inhomogeneous in an active star- 
forming galaxy (Sanders et al. 1988; Gordon, Calzetti & 
Witt 1997). In such an inhomogeneous ISM, optically-thin 
surface layers of optically-thick clumpy clouds are stripped 
by the radiation drag, and the stripped gas loses angular 
momentum, eventually accreting on to the center (Sato 
et al. 2004). Kawakatu & Umemura (2002) have shown 
that the inhomogeneity of ISM plays an important role 
for the radiation drag to attain the maximal efficiency. 
Then, the final mass of MDO is proportional to the total 
radiation energy from stars. Taking the realistic chemi- 
cal evolution into consideration, the radiation drag model 
predicts a mass ratio as Men/^buige — 0.001 in a galactic 
bulge (Kawakatu, Umemura, & Mori 2003). The theoret- 
ical upper limit of BH-to-bulge mass ratio is determined 
by the energy conversion efficiency of nuclear fusion from 
hydrogen to helium, i.e., 0.007 (Umemura 2001). Here, 
the question is whether the radiation drag mechanism can 
work also in small spheroidal systems like globular clus- 
ters. 

In this paper, we examine the possibility of the radi- 
ation drag-induced formation of massive BHs in globu- 
lar clusters. Here, we assume that the formation of a 
globular cluster begins with the coeval starburst within 
it. The outline of the physical processes we consider is 
shown in Figure 1: (I) First, the energy input from type 
II supernovae (SNe) drives the outflow of ISM from the 
system. Therefore, the system becomes optically-thin in 
an early evolutionary phase. (II) Afterward, intermediate 
mass (2 — 8M0) stars shed the gas envelope into the inter- 
stellar space, accumulating the ISM, and thus the optical 
depth of the system increases. In this phase, the radiation 



drag can work effectively and the mass accretion on to the 
center is induced. (HI) Finally, type la SNe, with type II 
SNe in newly formed stars, expel the ISM from the system 
again. Resultantly, the efficiency of the radiation drag de- 
scends abruptly. Through these physical processes, it is 
expected that the BH formation in a globular cluster is 
strongly related to the star formation history. 

In this paper, we attempt to elucidate whether the for- 
mation of massive BHs can be prompted by a radiation hy- 
drodynamic mechanism in globular clusters. Also, based 
on the present model, we reconsider the formation of mas- 
sive BHs in M15, Gl, uo Gen, the M33 nucleus, and the 
compact X-ray sources in M82. The paper is organized 
as follows. In Section 2, we build up the model for the 
chemical evolution and the radiation-hydrodynamic pro- 
cess. In Section 3, we investigate the relation between the 
star formation history and the final BH mass in a cluster. 
Then, we derive the condition on which globular clusters 
can possess massive BHs. In Section 4, we compare our 
predictions with the observational results for M15, Gl, uo 
Gen, the M33 nucleus, and the compact X-ray sources in 
M82, and discuss whether massive BHs can form through 
radiation-hydrodynamic process in these globular clusters. 
In Section 5, based on the present regime, we discuss ob- 
servational indications of the formation process of massive 
BHs in globular clusters. Section 6 is devoted to the con- 
clusions. 



(I) ISM is expelled due to (II) ISM ejected from intermediate (III) ISM is expelled due to 
supernovae. mass stars can be retained. supernovae again. 




optically thin optically thick optically thin 

Radiation drag is effective. 
So, the BH formation is promoted. 



Fig. 1. — The outline of the present scenario for the formation 
of a massive BH in a globular cluster. 

2. MODELS 

2.1. Chemical evolution of globular cluster 

We suppose a two-component system that consists of 
a spheroidal stellar cluster and dusty interstellar matter 
within it. First, we construct the model for the chemi- 
cal evolution of globular clusters, by using an evolution- 
ary spectral synthesis code 'PEGASE' (Fioc & Rocca- 
Volmerange 1997). At the initial epoch, we set a spherical 
system composed of gas with mass given in the range of 
Mo = 1O^~^M0. The radius, rg, is determined to follow a 
mass-radius relation, (rg/lOpc) = (Mq/IO^M©)^/^, found 
by Chiosi & Carraro (2002). It is assumed that the star 
formation starts with an overall starburst with a duration 
of ^ lO^yr, which reproduces the color-magnitude rela- 
tion of globular clusters at the present epoch (Yoshii & 
Arimoto 1987). We employ a stellar initial mass function 
(IMF) in the form of ^ = dn/d\ogm^ = ^(m^/M©)"", 
where is the stellar mass. The lower mass in the IMF 
is assumed to be m\ = 0.1 Mq and the upper mass to be 
rriu = 60 Mq. The index a is changed from 0.2 to 2.0 
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according to the previous analysis (Kawakatu & Umemura 
2002). The internal density in a cloud is assumed to be 
uniform, and the optical depth is determined by a dust-to- 
gas ratio that is calculated by TEGASE'. Thus, the total 
optical depth of globular cluster, tt , is a function of time. 
Also, the optical depth of a cloud and therefore the overall 
optical depth of globular cluster depend on the cloud size 
Tc- Here, = O.Olrg is assumed as a fiducial case. But, 
we have confirmed that no essential difference in the final 
BH mass is found by changing Tc so as to enhance tt by 
an order. 



2.2. Mass accretion due to radation drag 



Next 
cluster. 



C{t) 



we model the formation of a BH in a globular 
based on the radiation drag-driven mass accre- 
tion. The radiation drag, which drives the mass accretion, 
originates in the relativistic effect in absorption and subse- 
quent re-emission of the radiation. This effect is naturally 
involved in relativistic radiation hydrodynamic equations 
(Umemura, Fukue, & Mineshige 1997; Fukue, Umemura, 
& Mineshige 1997). The angular momentum transfer in 
radiation hydrodynamics is given by the azimuthal equa- 
tion of motion in cylindrical coordinates. 



1 d{rv^) _ Xd 



dt 



Ad 



{E + P^^)i 



(2) 



in order to see the effects of the slope of IMF on the final 
results. Also, stars with masses greater than SMq are pos- 
tulated to undergo type II supernova explosions. As for 
the type la supernova, an evolutionary spectral synthe- 
seis code TEGASE' adopts the single-degenarate scenario 
(Nomoto et al. 1994) for the type la supernova progenitor, 
where a white dwarf (WD) in a close binary undergoes a 
thermonuclear explosion when the companion star evolves 
off the main sequence and transfers a large enough amount 
of mass to the WD. ^ From the evolution model of pro- 
genitor stars, Hachisu et al. (1996, 1999) obtained the 
mass range of companion stars of the WDs that become 
type la supernovae as [O.9M0, S.SM©]. The star formation 
rate (SFR) per unit mass , C(t), is assumed to be in pro- 
portion to the gaseous mass fraction, fg{t) = Mg(t)/Mo, 
where Mg{t) is the total gas mass at time t. The SFR is 
given by 

^ Pfg t < twl , 

kfg (twl < t < tw2), (1) 
10 (t>tw2). 

The epoch t^i is determined by the age of the minimal 
star (SMq) that undergoes a type II supernova, that is, 
^wi = 6 X lO^yr. Before twi, the interstellar gas is assumed 
to blow out from the system due to the energy input from 
successive type II SNe. tw2 is the second blowout epoch 
when the interstellar gas is expelled by type la SNe as well 
as type II SNe in newly formed stars. The second blowout 
epoch is estimated by the epoch when the kinetic energy 
(£'kin) equals the binding energy (E'bin)- The kinetic en- 
ergy is evaluated by E'kin = elO^^A^sNerg, where A^sn is 
the number of SNe, which is calculated by TEGASE'. e 
is the efficiency at which the total SN energy is converted 
to the kinetic energy of blast wave. The Sedov solution 
gives e = 0.28 for the adiabatic exponent of 7 = 5/3. Af- 
ter tw2, the system is assumed to be ISM-deficient and 
thus optically thin. The coefficient P is assumed to be 
P ^ (10''yi")~^(Yoshii & Arimoto 1987), and k is changed 
in order to see the dependence on the SFR. 

As for the stellar component, we assume star clusters 
with a specific stellar initial mass function (see above). 
Stars are distributed uniformly inside the spheroidal sys- 
tem. Actually, A^*(= 300) small star clusters are dis- 
tributed randomly. In addtion, the angular momentum 
is smeared, assuming rigid rotation, according to the spin 
parameter of A = {Jt\Et\^/^)/{GM^^'^) = 0.05, where Jt, 
E't, and Mt are respectively the total angular momentum, 
energy, and mass (Barnes & Efstathiou 1987; Heavens & 
Peacock 1988; Steinmetz & Bartelmann 1995). 

As for ISM, we consider clumpy matter, since the ISM 
is observed to be highly inhomogeneous in active star- 
forming regions (Sanders et al. 1988; Gordon, Calzetti, 
& Witt 1997). Here, Nc{= 10"^) identical clouds are dis- 
tributed randomly. (It is noted that simulations with a 
three times larger number of clouds did not lead to any 
fundamental difference in the final BH mass, although at 
least 10^ clouds are necessary to treat the radiative trans- 
fer effect properly in clumpy ISM.) The distribution and 
motion of ISM is the same as those of the stellar compo- 
nent. We assume that the cloud covering factor is unity 

^ In present-day galaxies, it is well known that more than half of main-sequence stars are observed in binaries, while it is not observationally 
clear if Population III stars are able to form in binary systems. But, Saigo, Matumoto & Umemura (2004) have recently shown by numerical 
simulations that a significant fraction of Population III stars are expected to form in binary systems, and thus type la supernovae would occur 
even in the metal- deficient environments. 



where E is the radiation energy density, is the radiation 
flux, P^^ is the radiation stress tensor, and Xd is the mass 
extinction coefficient which is given by Xd = nd(Jd/pgas 
with the number density of dust grains na, the dust cross- 
section (Jd and the gas density Pgas- By solving radiative 
transfer with including dust opacity, we evaluate the ra- 
diative quantities, E^ F^, and P"^"^, and thereby obtain 
the total angular momentum loss rate. Then, we can es- 
timate the total mass of the dusty ISM accreted on to a 
central massive drak object, Mmdo- By using the rela- 
tion, Mgas/^gas = —J/Ji where J and Mgas are the total 
angular momentum and gas of ISM. Then, the MDO mass 
is assessed as 



Ml 



MDO 



= / Mg^dt = - I 
Jo Jo 



M gas -J dt. 



(3) 



In the optically-thick regime of the radiation drag, 
Mgs^sJ/J = J^[Lstar,y{t)/c^]diy, where Lstar,!. is the to- 
tal luminosity of globular cluster at frequency u. The ra- 
diation drag efficiency depends on the optical depth r in 
proportion to (1 — e~^) (Umemura 2001). Thus, the total 
mass of MDO can be expressed by 



Jo Jo 



■At) 



(1 - e-^-W) 



dudt, 



(4) 

where r^, is the optical depth of the globular cluster mea- 
sured from the center. Here, we estimate the evolution 
of by using TEGASE'. The efficiency 

^drag is found to 

be maximally 0.34 (Kawakatu & Umemura 2002). After 
tw2, the radiation drag is inefficient, because the system 
becomes ISM-deficient and optically thin. In this paper. 
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at t < twi and t > t^2, is assumed to drop abruptly to 
r^y <C 1. Hence, the final mass of an MDO is given by 

Mmdo = r?drag r ^^^^^ (l - e--(*)) dvdt. 

Jt^i Jo ^ ^ ^ 

(5) 

2.3. The fate of the MDOs 

In the previous section, we estimated the MDO mass in 
the context of the radiation drag-induced mass accretion. 
However, the MDO itself does not evolve into the mas- 
sive BH directly, because the radiation drag is not likely 
to remove the angular momentum thoroughly (Sato et al. 
2004). Hence, some residual angular momentum will ter- 
minate the radial contraction of the accreted gas. Thus, 
the MDO is likely to be a compact rotating disk. We 
should consider the further collapse of the MDO through 
other physical mechanisms. In the MDO, the viscosity 
is expected to work effectively, because the timescale for 
viscous accretion is shortened by the radiation drag (Mi- 
neshige, Tsuribe, & Umemura 1998). Thus, the MDO 
would be a massive self-gravitating viscous disk. For a 
massive self-gravitating viscous disk, some self-similar so- 
lutions are known to give an inside-out disk collapse in 
a fiat temperature distribution of a disk (Mineshige & 
Umemura 1996, 1997; Tsuribe 1999). A fiat temperature 
profile is plausible as long as the viscous heating is bal- 
anced with the radiative cooling in the geometrically thin 
disk (see Mineshige & Umemura 1996 for details). In fact, 
Tsuribe (1999) derived a series of self-similar solutions for 
rotating isothermal disk, taking into account the growth 
of the central point mass, and then provided a convenient 
formula for the inside-out mass accretion rate for a variety 
of outer flows. 



where a is the viscous parameter, Cg is the sound speed 
and Q is theToomre's Q which is kc^/tiGT, for the epicy- 
cle frequency and the surface density E. Tsuribe (1999) 
has found that an accretion with Q ?^ 2 is stable for a wide 
range of a. The critical accretion rate is given by 

M„^|=0.24Meyr-i(^)'^'. (7) 

It must be noted that the viscous accretion rate is sen- 
sitively dependent on the disk temperature. Under an 
intense starburst, the disk is exposed to the strong ul- 
traviolet radiation. Then, the disk could be heated up to 
lO^K, although the detail is dependent on the dust ex- 
tinction and radiative cooling. Here, we employ lO^K as 
the disk temperature. Through this inside-out collapse, 
almost all the MDO can accrete on to a centeral core. Fi- 
nally, the core is likely to be a rigidly-rotating very massive 
star (VMS), because the angular momentum transfer via 
viscosity works to smear out any differential rotation in a 
self-gravitating system. Such VMSs pessess large helium 
cores that reach carbon ignition. The fate of VMSs sensi- 
tively depends on the initial mass (Bond, Arnett, & Carr 
1984; Heger & Woosley 2002). A VMS below 26OM0 (in 
the range 140 — 26OM0) results in a pair-instability super- 
nova, leaving no remnant, because the core of such a star 
undergoes the electron-positron pair creation instability 



after the helium burning (Barkat, Rakavy, & Sack 1967; 
Bond, Amett, & Carr 1984; Kippenhahn & Weigert 1990; 
Flyer, Woosely, & Herger 2001). Above - 26OM0, a VMS 
collapses directly into a BH because the nuclear burning 
cannot halt the core collapse (Fryer et al. 2001; Heger et 
al. 2003). There must be the effect of rotation, because 
the progenitor stars are likely to be rotating (Bond et al. 
1984; Fryer et al. 2001). Recently, Shibata (2004) finds 
that a rigidly-rotating VMS with several IOOMq can be 
unstable in general relativity for softer equations of state, 
eventually forming a BH. 

Here, we should pay attention to the metallicity depen- 
dence of the VMS formation. In the solar abundance, any 
stars more massive than lOOM© are not expected to form 
(Fryer et al. 2001; Heger et al. 2003). On the other 
hand, VMSs can preferentially form in metal- free gas (e.g 
Nakamura & Umemura 2001, and references therein). The 
problem is the critical metallicity, under which VMSs can 
form. Recently, Omukai & Palla (2003) have explored this 
problem in detail, and found the critical metallicity to be 
around O.OIZ©, under which dust opacity is not effective to 
hinder the mass accretion onto a protostar. Hence, in the 
light of the critical metallicity, VMSs are expected to form 
in globular clusters. Taking such theoretical advances into 
account, we assume 26OM0 as the minimum mass for the 
formation of massive BHs. If Mmdo exceeds this thresh- 
old mass (26OM0), a massive BH whose mass is equal to 
a MDO mass would form in a globular cluster. However, 
it is not clear enough whether the self-similar solutions 
employed in this paper are feasible in realistic situations 
for the MDO formed via the radiation drag. The problem 
should be explored with sophisticated numerical simula- 
tions, which are left for the future study. In this paper, 
we regard a BH mass estimated here as a maximal mass 
of BH within a globular cluster. 

3. RESULTS 

By changing physical parameters a and /c, which control 
the star formation, we investigate the relation between the 
formation of massive BHs and the star formation history. 
First, we examine the dependence on the stellar IMF and 
then on the star formation rate. Finally, we derive the 
relation between the mass of globular cluster and the BH 
mass. 

3.1. Dependence on Initial Mass Function 

Here, we investigate the effects of the IMF slope {a) 
on the final BH mass, with settling Mq = lO^M© and 
k = 8.6Gyr"^ In Fi gure 2, the resultant mass of BH 
(Mbh) is shown for the range ofa = 0.2 — 2.0. The details 
of the results are summarized in Table 1 . As seen in Figure 
2, we find that the mass of central massive objects exceeds 
the threshold mass (26OM0) for the BH formation in the 
range of a = 0.2 — 2.0. Also, the BH mass shows a peak 
around c± 1. The reason can be understood as follows. 
First, let us remind that the final BH mass is determined 
by the optical depth of the system and the duration be- 
tween twi and tw2, as shown in equation (5). In the present 
scenario, the optical depth is controlled by the mass loss 
from stars with < SMq. If a is smaller, a larger number of 
massive stars form and therefore the optical depth can be 
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increased in a shorter timescale by the mass loss from in- 
termediate mass stars. But, simultaneously an increasing 
number of SN explosions shorten the second blowout time 
tw2 • Table 1 shows that reduces vastly with decreasing 
a, so that Mbh lessens. On the other hand, if a is larger, 
the fraction of massive stars decreases. Result antly, fewer 
SN explosions lengthen t^2- But, the optical depth can- 
not augment owing to a reduced number of intermediate 
mass stars. As a consequence, Mbh lessens with increas- 
ing a. By these two effects, a 1 leads to the maximum 
efficiency of the radiation drag. 

The change of a brings a different mass-to-luminosity 
ratio and metallicity of a globular cluster. The shaded 
region satisfies the observed mass-to-luminosity ratio of 
globular clusters, 1.0-3.0 in solar units (e.g., Mandushev 
et al. 1991; Djorgovski et al. 2003), and also the metal- 
licity, ^ 0.01 - 0.1^0 (Lee 1990; Chaboyer, Demarque & 
Sarajedini 1996), where Zq is the solar metallicity. If the 
slope of IMF is in the range of a = 0.7 — 1.7, the properties 
of globular clusters are consistent with the observations. 
In this range of ce, massive BHs with ^ 10^ Mq are ex- 
pected to form. 

In this calculation, we have assumed m\ = O.IM© and 
TTiu = 60Mq. To see the effect of a different choice of the 
lower and upper mass of the IMF, we also investigated 
the case of (mi, mu)=(O.5M0, GOM©), (O.IM®, 3OM0), 
and (O.IM0, 12OM0), all of which reproduce the observed 
mass-luminosity ratio. As a result, we found that the fi- 
nal BH mass is not strongly dependent on m\ and rriu. 
To conclude, as long as the mass-to-luminosity ratio and 
metallicity satisfy those observed, the final BH mass is not 
sensitive to the stellar initial mass function. 
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Fig. 2. — The dependence of the final BH mass (Mbh in units of 
Mq) on the IMF slope (a). Here, Mq = lO^M© and k = 8.6Gyr-^ 
are assumed. The solid line shows our prediction. The shaded region 
satisfies to the observed mass-to-luminosity ratio and metallicity in 
globular clusters. The BH mass is maximal around o; ~ 1. 

3.2. Dependence on Star Formation Rate 

In this section, to examine the dependence of the BH 
mass on the star formation history in globular clusters, we 
alter the coefficient k from O.lGy"^ to 300Gy~^. Here, we 
assume Mq = lO^M©, and set a = 0.95 which gives the 
maximal radiation drag efficiency for the BH formation, 
as shown in the previous section. In Figure 3, the final 
BH mass is shown against k. The corresponding star for- 
mation timescale, which is defined by tsF = is also 
shown on the upper horizontal axis. The maximal optical 



depth at U -band is attached on each point simulated. The 
shaded region satisfies the observed mass-to-luminosity ra- 
tio and metallicity of globular clusters. As seen in Figure 
3, if tsF is shorter than ~ lO^yr, the BH mass becomes 
lower. This is because tsF is shorter than the typical age 
of intermediate mass stars, ^ lO^yr, the optical depth 
cannot increase to a high level before t^2- In contrast, 
if tsF is longer than t*, massive BHs with > lO^M© can 
form, because the optical depth of the system can reach 
a high level. However, if tsF is ^ lO^^yr, the number of 
intermediate mass stars increase slowly with time. Hence, 
the BH cannot grow rapidly. Hence, the BH mass is al- 
most saturated for longer tsF- As a result, we find that, as 
long as the mass-to-luminosity ratio and metallicity satisfy 
those observed, the final BH mass is almost regardless of 
the star formation rate, and massive BHs with ^ IO^Mq 
are expected to form. 
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Fig. 3. — The final BH mass (Mbh) against the star formation 
rate coefficient (k), assuming Mq = 10^ Mq and a = 0.95. The 
abscissa is the coefficient in units of Gyr"-*^, and the ordinate is the 
ffnal BH mass in units of solar mass. The upper horizontal axis is 
the star formation timescale, tsF in units of yr. The attached values 
are the maximal optical depth at [/-band. For tsF > lO^yr, the BH 
mass is almost saturated. For tsF < lO^yr, the BH mass decrease 
rapidly, because the system becomes optically thin. The shaded re- 
gion satisffes the observed mass-to-luminosity ratio and metallicity 
in globular clusters (Z* = 0.01 - O.IZq). 

3.3. Globular Cluster-to-BH Mass Relation 

Finally, we derive the relation between the final BH 
mass, Mbh, and the final stellar mass of globular clusters, 
Mstar, final- Here, we set k to be constant as = 8.6Gyr~^C 
but a is changed. We consider the range of initial mass as 
Mo = 1O^~^M0, because we focus on the final mass range 
of Mstar, final = lO^^^M©. The rcsults are shown in Fig- 
ure 4. As shown in previous sections, the final BH mass 
is almost independent of the stellar initial mass function 
and star formation rate, as long as the resultant proper- 
ties of globular clusters are consistent with the observa- 
tions. In Figure 4, the shaded area is the region which 
satisfies the observed mass-to-luminosity ratio and metal- 
licity of globular clusters. The region of no BH means 
that the mass of MDO is less than 26OM0, which is the 
present criterion for a massive BH. Figure 4 shows that 
massive BHs can form only for the final mass of glob- 
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ular clusters with > 6 x lO^M©. The putative BH-to- 
bulge mass relation, Mbh = O.OOlMbuige, is shown by a 
dashed line. We can see that the predicted BH-to-globular 
cluster mass ratio (Mbh /^star, final) is considerably lower 
than the BH-to-bulge mass ratio. The main reason is that 
the duration when the system is optically thick is shorter 
in globular clusters than that in galactic bulges, because 
the supernova explosions are more devastating for smaller 
spheroidal systems like globular clusters. It is also noted 
that the relation between the BH mass and the globular 
cluster mass is not linear, in contrast to the BH-to-bulge 
mass relation. 




105 106 107 108 

Mstar/mal [^^O 1 

Fig. 4. — The relation between the final BH mass (Mbh) and 
the final stellar mass of globular cluster (Mgtar final)- Here, we as- 
sume k = S.eGyr"-*^. The ordinate is the BH mass in units of solar 
mass, and the abscissa is the final mass of globular cluster in units 
of solar mass. The thick solid lines show the results for diflFerent 
initial masses of the system, Mq. The attached values denote the 
IMF index a. The shaded area shows the region which satisfies the 
observed mass-to-luminosity ratio and metallicity. The BH-bulge 
mass relation, Mbh = O.OOlMbuige? is shown by a dashed line. The 
predicted BH fraction is considerably smaller than 10~^. 

4. COMPARISON WITH OBSERVATIONS 

To compare the present prediction with the observations 
on the massive BHs in globular clusters, we translate the 
obtained Mbh — Mgtar, final relation into Mbh — cr relation, 
based on the virial theorem, where a is the stellar velocity 
dispersion in globular clusters. The resultant Mbh — cr 
relation is shown in Figure 5. Here, we compare the pre- 
diction with BH candidates in Ml 5, Gl, oj Gen, the M33 
nucleus, and the compact X-ray sources in M82. 

In a globular cluster Ml 5 in the Milky Way, the BH 
mass is estimated to be Mbh = i^-^tl r) x lO^M© by 
Gerssen et al. (2003), which includes the possibility of 
no BH. By the detailed comparison with numerical sim- 
ulations, Baumgardt et al. (2003a) claim that a massive 
BH with > 500 — lOOOM© is not indispensable to account 
for the observational data in Ml 5. The total stellar mass 
of M15 is roughly IO^Mq and the velocity dispersion is 
a ~ 14km/s. As seen in Figure 5, the present model pre- 
dicts that no massive BH forms for the velocity dispersion 
in M15. Hence, the prediction of the present model sup- 
ports the case of no BH. 

In a globular cluster Gl in M31, a candidate for a mas- 
sive BH with Mbh 2 x 10"^ M© is reported (Gebhardt, 
Rich & Ho 2002). But, Baumgardt et al. (2003b) argue 
that the observational data in Gl can be well fitted even 
without a massive BH. Hence, the above mass of BH in Gl 



should be considered as an upper limit. The total mass of 
Gl is (0.7 — 2) X 10^ Mq and the velocity dispersion is esti- 
mated as (7 ~ 25km/s (Meylan et al. 2001). It is predicted 
for this velocity dispersion that no massive BH forms if the 
system mass is initially lower than 5 x lO^M©, while a BH 
with ^ 5OOM0 can form if the system mass is initially 
higher than that mass. 

In addition, uj Gen possibly harbors a massive BH. uj 
Gen is a globular cluster in the Milky Way in which mem- 
ber stars show a wide range of metallicity. The total stellar 
mass of UJ Gen is 5 x lO^M©, and the velocity dispersion 
is cr ?^ 22km/s (Meylan et al. 1995). As seen in Figure 5, 
whether co Gen can have a massive BH with ^ 3OOM0 is 
marginal. Latest radio observations have suggested that 
the black hole mass in co Gen should be less than about 
IOOM0 for the spherical Bondi-Hoyle accretion rate, but 
the data are marginally consistent with a black hole of 
about IOOOM0 for more plausible accretion rate (Mac- 
carone et al. 2005). 
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Fig. 5. — Comparison of the theoretical prediction with the ob- 
servational data on BH mass in globular clusters. The ordinate is 
the BH mass (Mbh) in units of solar mass, and the abscissa is the 
velocity dispersion of globular cluster (cr) in units of km/s. The 
parameters are the same as Figure 4. The shaded region is the pre- 
diction of the present analysis. The relation between the BH mass 
and the velocity dispersion of galactic bulges is shown by a dotted 
line, Mbh = 1-3 x loSM0(cr/2OOkm/s)^-O2 (Tremaine et al. 2002). 
The observational data of M15, Gl, ix) Cen, M33 nucleus and M82 
are plotted. 

On the other hand, no strong evidence of a massive BH 
is found in a stellar cluster at the center of M33 (Merritt, 
Ferrarese, & Joseph 2001; Gebhardt et al. 2001). The M33 
nucleus has the velocity dispersion of cr = 21 — 34km/s 
(Merritt, Ferrarese, & Joseph 2001). For this velocity 
dispersion, our model predicts that a massive BH with 

lO^M© can form only if the initial system mass is higher 
than 5 x lO^M©. 

Finally, Matsumoto et al. (2001) have identified nine 
bright compact X-ray sources in M82 by Chandra obser- 
vatory. The brightest source (No. 7 in their Table 1) has 
a luminosity of 9 x lO^^ergs"^. If this luminosity is as- 
sumed to be the Eddington luminosity, the derived mass 
of a black hole is 35OM0. From the infrared luminosity, the 
total mass of the cluster is estimated to be ~ 5 x lO^M© 
and the velocity dispersion is to be llkm/s (McCrady et 
al. 2003). As seen in Figure 5, the present model predicts 
that no BH forms for this velocity dispersion, although the 
possibility of BH formation via runaway stellar collision is 
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proposed (Portegies Zwart & McMillan 2002; Portegies 
Zwart et al. 2004). 

5. DISCUSSION 

In the previous sections (§3 and §4), we have shown 
that massive BHs are expected to be hosted by globular 
clusters. Here, we consider the detect ability of the BH 
formation in globular clusters. The present scenario of 
BH formation is summarized as follows: (i) First, a mas- 
sive self-gravitating viscous disk (MDO) forms owing to 
the radiation drag in a timescale of ~ 10^ — lO^yr, and 
the core evolves into a VMS. (ii) Secondly, the VMS col- 
lapses directly into a BH if its mass is higher than 26OM0. 
(iii) Finally, the BH grows to a massive BH by the mass 
accretion from the ambient gas. 

A VMS more massive than 26OM0 evolves in a short 
lifetime of < lO^yr (e.g., Schaerer 2002). Therefore, the 
stages (i) and (ii) correspond to a high redshift, if we 
consider the age of globular clusters. Hence, it would be 
hard to detect a VMS itself, because the apparent magni- 
tude is too faint to observe. However, a VMS may result 
in a gamma-ray burst (GRB), if it is Population H and 
a forming BH entails the sufficient angular momentum 
(Heger et al. 2003). GRBs are detectable even at red- 
shift higher than 10 (Lamb & Reichart 2001). The present 
scenario predicts that GRB events would take place at 
the centers of forming globular clusters at early stages of 
^ 10^ - lOVr- 

In the stage (iii), the luminosity of BH accretion is ex- 
pected to be L = TjMc^^ where M is the mass accretion 
rate and 77 is the efficiency. The temperature in the inner- 
most regions of accretion disk and the efficiency r] depend 
on whether M is higher or lower than the Eddington ac- 
cretion rate Me = L^/c^ (e.g. Kato, Mineshige, & Fukue, 
1998), where Le is the Eddington luminosity given by 

L..L25xlO"e,..-'(^). (8) 

If M ^ Me, a so-called standard disk forms. Then, the 
efficiency is 77 0.1 and the temperature is given by 

where rgch is the Schwarzschild radius. If M ^ Me, the so- 
lution of accretion is optically thin and called a radiatively 
inefficient accretion flow {RIAF). Then, r] ^ 0.1(M/Me) 
and the temperature is on the order of 10^ K for electrons. 
If the accretion rate is super-Eddington (M > Me), the 
solution is optically-thick and called a slim disk, where the 
temperature is higher than the standard disk. In this type 
of disk, the photon-trapping effect in the accretion flow 
plays an important role and then the efficiency becomes 
T] ^ 0.1(M/Me)~^/^. Actually, the accretion luminosity 
can achieve up to ^ 7Le (Ohsuga et al. 2002, 2003). 

Thus, if M > Me, it is expected that an accreting BH in 
a globular cluster would be a far-UV and X-ray source with 
L > I/E- If a slim disk is considered, L « lO^^ergs"^ is 
realized even for a BH with several lOOM©. Such a source 



may be observed as a ultraluminous X-ray source (ULX), 
the luminosity of which is 10^^ — lO^^ergs"^. The lifetime 
of the accretion phase is pretty shorter than the present 
age of globular clusters. But, in some globular clusters, 
the final accretion phase may be detected. In fact, Liu & 
Bregman (2005) have recently reported that some ULXs 
are located at old globular clusters. These ULXs can be 
the candidates for accreting BHs with lO^M© formed by 
the radiation hydrodynamic process. 

6. CONCLUSIONS 

Based on the radiation hydrodynamic model, we have 
examined the possibility for the formation of massive BHs 
in globular clusters, with incorporating the realistic chem- 
ical evolution. In the present model, the mass of a central 
massive dark object in a globular cluster is determined by 
the duration of the optically- thick phase. We have found 
that massive BHs can form only in globular clusters with 
the the final mass of > 6 x lO^M© or the velocity disper- 
sion of > 20km/s. The BH mass is almost independent of 
the stellar initial mass function and star formation rate. 
The results show that massive BHs are not likely to form 
in typical globular clusters with ^ IO^Mq. The predicted 
BH-to-globular cluster mass ratio (MBn/Mstar, final) is con- 
siderably lower than the BH-to-bulge mass ratio. 

We have applied the present model to some globular 
clusters, in which the possibilities of massive BHs are re- 
ported. In M15, no massive BH is predicted to form, 
although the evidence of a massive BH is under debate. 
Gl can harbor a massive BH with ^ SOOMq, which 
is smaller by an order of magnitude than an estimate, 
Mbh ^ 2 X lO^Mo, by Gebhardt, Rich & Ho (2002). In 
Lj Gen, the formation of a massive BH with ^ 3OOM0 is 
marginal. As for M82, it is suggested that massive BHs 
with 35OM0 can reside in bright compact X-ray sources. 
But, no BH is expected to form in the present radiation 
hydrodynamic model. Finally, the present scenario pre- 
dicts that a final phase of BH growth in a globular cluster 
can be observed as a ULX with ~ lO^^ergs"^. Also, GRBs 
may arise in early phases (10^ — lO^yr) of globular cluster 
formation. 

In this paper, we have investigated the possibility of the 
BH formation from a hydrodynamic point of view. There 
may be a different path to form massive BHs in globu- 
lar clusters. To eclucidate the mechanism for the forma- 
tion of massive BHs, the high resolution observations due 
to the mult i- wavelength, especially optical, far-UV, X-ray 
and gamma-ray on the central regions of globular clusters 
would be greatly important in the future. 
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Table 1. BH mass-to-IMF Slope in Globular Cluster 
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(a) Mgtar, final is the mass of globular cluster in the present day. (b) ru,max represents the maximal optical depth at C/-band. 



REFERENCES 



Barkat, Z., Rakavy, G., Sack, N. 1967, Phys. Rev. Lett., 18, 379 

Barnes, J., & Efstathiou, G. 1987, ApJ, 319, 575 

Baumgardt, H., Hut, P., Makino, J., McMillan, S., & Portegies 

Zwart, S. 2003a, ApJ, 582, L21 
Baumgardt, H., Hut, P., Makino, J., McMillan, S., & Portegies 

Zwart, S. 2003b, ApJ, 589, L25 
Bond, J. R., Arnett, W. D., &; Carr, B. J. 1984, ApJ, 280, 825 
Chaboyer, B., Demarque, P., & Sarajedini, A. 1996, ApJ, 459, 558 
Chiosi, C., &; Carraro, G., MNRAS, 2002, 335, 335 
Djorgovski, S.G.et al., "New Horizons in Globular Cluster 

Astronomy", eds. G. Piotto, G. Meylan, S.G. Djorgovski, and M. 

Riello, 2003, ASPCS, in press 
Fioc, M., &Rocca-Volmerrage, B. 1997, A&A, 326, 950 
Fryer, C. L., Woosley, S. E., & Heger, A. 2001, ApJ, 550, 372 
Fukue, J., Umemura, M., & Mineshige, S. 1997, PAS J, 49,673 
Gebhardt, K., Rich, R.M., &; Ho, L. 2002, ApJ, 578, L41 
Gebhardt, K. et al. 2001, AJ, 122, 2469 
Gerssen, L. et al. 2003, AJ, 125, 376 

Gordon, K., Calzetti, D., & Witt, A.N. 1997, ApJ, 487, 625 
Hachisu, I., Kato, M., & Nomoto, K. 1996, ApJ, 470, L97 
Hachisu, L, Kato, M., & Nomoto, K. 1999, ApJ, 522, 487 
Heaveans, A., & Peacock, J. 1988, MNRAS, 232, 339 
Heger, A., Fryer, C.L., Woosley, S.E., Langer, N., & Hartmann, 

D.H., 2003, ApJ, 591, 288 
Heger, A., &; Woosley, S.E., 2002, ApJ, 567, 532 
Hut, P., et al. 1992, PASP, 104, 981 

Kato, S., Mineshige, S., &; Fukue, J. 1998, Black Hole Accretion 

Disks (Kyoto: Kyoto Univ. Press) 
Kawakatu, N., Umemura, M., & Mori, M. 2003, ApJ, 583, 85 
Kawakatu, N., & Umemura, M. 2002, MNRAS, 329, 572 
Kippenhahn, R., & Weigert, A. 1990, Stellar Structure and 

Evolution (Berlin: Springer) 
Kormendy, J., &; Richstone, D. 1995, ARA&A, 33, 581 
Lamb, D.Q., &; Reichart, D.E. 2001, ApJ, 536, 1 
Lee, H.M., 1987, ApJ, 319, 801 
Lee.Y.-W., 1990, ApJ, 363, 159 
Liu, J., & Bregman, J. 2005, ApJ, 157, 59 



Maccarone, T. J., Fender, R. P., &; Tzioumis, A. K. 2005, MNRAS, 
356, L17 

Mandushev, G., Spassova, N., &; Staneva, A. 1991, A&A, 252, 94 
McCrady, N., Gilbert, A.M., Graham, J.R., 2003, ApJ, 596, 240 
Matsumoto, H. et al. 2001, ApJ, 547, L25 

Merritt, D., Ferrarese, L., & Joseph, C. 2001, Science, 293, 1116 
Meylan, G. et al. 2001, AJ, 122, 830 
Meylan, G. et al. 1995, A&A, 303, 761 

Mineshige, S., Tsuribe, T., &; Umemura, M. 1998, PAS J, 50, 233 

Mineshige, S., & Umemura, M. 1996, ApJ, 469, L49 

Mineshige, S., & Umemura, M. 1997, ApJ, 480, 167 

Nakamura, F. & Umemura, M. 2001, ApJ, 548, 19 

Nomoto, K., Yamaoka, H., Shigeyama, T., Kumagai, S., & 

Tsujimoto, T. 1994, in Les Houches Session LIV, ed. S.A. Bludman 

(Amsterdam: Elsevier), 199 
Ohsuga, K., Mineshige, S., Mori, M., & Umemura, M. 2002, ApJ, 

574, 315 

Ohsuga, K., Mineshige, S., &; Watarai, K. 2003, ApJ, 596, 429 
Omukai, K., & Palla, F. 2003, ApJ, 589, 677 

Portegies Zwart S.F., Baumgardt, H., Hut, P., Makino, J., & 

Mcmihan, S.L.W., 2004, Nature, 428, 724 
Saigo, K., Matsumoto, T., & Umemura, M., 2004, ApJ, 615, L65 
Sanders D.B., Soifer B.T., Elias J.H., Madore B.F., Matthews K., 

Neugebauer G., Scoville N.Z., 1988, ApJS, 325, 74 
Sato, J., Umemura, M., Sawada, K., & Matsuyama, S. 2004, 

MNRAS, 354, 176 
Schaerer, D. 2003, A&A, 382, 28 

Shakura, N. I., & Sunyaev, R. A.. 1973, A&A, 24, 337 

Shibata, M. 2004, ApJ, 605, 350 

Steinmetz, M., & Bartelmann, M. 1995, 272, 570 

Tsuribe, T., & Umemura, M. 1997, ApJ, 486, 48 

Tsuribe, T. 1999, ApJ, 527, 102 

Tremaine, S., et al. 2002, ApJ, 574, 740 

Quinlan, G.D., &; Shapiro, S.L. 1987, ApJ, 321, 199 

Umemura, M. 2001, ApJ, 560, L29 

Umemura, M., Fukue, J., &; Mineshige, S. 1997, ApJ, 479, L97 
Umemura, M, Loeb, A., & Turner, E. L. 1993, 419, 459 
Yoshii, Y., & Arimoto, N. 1987, A&A, 188, 13 



